Theory of Phonon Hall Effect in Paramagnetic Dielectrics 
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Based upon spin-lattice interaction, we propose a theoretical model for the phonon Hall effect in 
paramagnetic dielectrics. The thermal Hall conductivity is calculated by using the Kubo formula. 
Our theory reproduces the essential experimental features of the phonon Hall effect discovered 
recently in ionic dielectric Tb3GasOi2, including the sign, magnitude and linear magnetic field 
dependence of the thermal Hall conductivity. 
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When an electrical current flows through a conduc- 
tor with direction perpendicular to an applied magnetic 
field, a transverse electrical current may be generated in 
the third perpendicular direction. This is well known as 
the Hall effect, and is due to the electromagnetic Lorentz 
force on the charge carriers. Accompanied transverse 
heat current also flows in the conductor, simply because 
the charge carriers carry energy, known as the Righi- 
Leduc effect. Two interesting variants of the conven- 
tional Hall effect is the anomalous Hall effect (AHE) in 
ferromagnets |l| and the spin Hall effect in nonmagnetic 
conductors 0, 0] , where the electron spin-orbit coupling- 
plays an essential role. The AHE is characterized by 
an anomalous contribution to the Hall resistivity com- 
ing solely from the magnetization. Intuitively, one would 
not expect a Hall effect for phonons, which do not carry 
charges and do not couple to the magnetic field directly. 
Remarkably, by applying a magnetic field perpendicular 
to a heat current flowing through a sample of the param- 
agnetic material TbsGasO^, Strohm, Rikken and Wyder 
observed very recently a temperature difference of up to 
200 microkelvin between the sample edges in the third 
perpendicular direction |^|. Since TbaGasO]^ is a dielec- 
tric, and so the Righi-Leduc effect can be ruled out. The 
temperature difference is attributed to the phonon Hall 
effect (PHE) jj], which becomes another intriguing and 
puzzling phenomenon in solid state physics. 

At the experimental low temperature 5.45K 0, exci- 
tation of optical phonons is unlikely, and thermal con- 
duction should be carried by acoustic phonons. While 
Tb3Ga50i2 is an ionic material, in a perfect lattice, each 
unit cell must be charge neutral. In the acoustic phonon 
modes, each unit cell vibrates as a rigid object without 
relative displacements between its constituting atoms , 
and hence does not acquire a net Lorentz force in a mag- 
netic field. Theoretical understanding of the physical 
mechanism underlying the PHE is so far absent. 

In this Letter, we propose a theoretical model based 
upon the Raman spin-lattice interaction for the PHE. 
The PHE is discussed to be a phonon analogue to the 
AHE. The thermal Hall conductivity of the phonons in 



the clean limit is calculated by using the Kubo formula. 
The theory can explain the essential features of the exper- 
imental data for Tb3Ga50i2, including the sign, magni- 
tude and linear magnetic field dependence of the thermal 
Hall conductivity. 

We consider a sample of a paramagnetic dielectric with 
volume V, which has a cubic lattice structure with sym- 
metry axes parallel to the axes of the coordinate system. 
For simplicity, we assume that only one rare-earth ion 
in each unit cell is paramagnetic. The relevant model 
Hamiltonian is written as 



H 



(1) 



where a qr7 creates an acoustic phonon of wavevector q 
and polarization a. The well-known Debye model is used 
to describe the acoustic phonons. We designate u = 
for the longitudinal acoustic (LA) phonons with disper- 
sion relation uj q o = c^q, and a = 1 and 2 for the trans- 
verse acoustic (TA) phonons with u> q i — uj q 2 — c T q. 
V represents the interaction between the phonons and 
the electronic spins and orbital an gula r momenta of the 
paramagnetic ions 0, HUll, often called the 
spin-lattice interaction. For Tb3Ga50i2, the rare-earth 
ions Tb + are paramagnetic with large magnetic mo- 
ments and may be responsible for the spin-lattice in- 
teraction. Since microscopic calculation of the spin- 
lattice interaction is difficult, phenomenological descrip- 
tion based upo n sy mmetry considerations is usually em- 
ployed 0)00) El EI- Two most important spin-lattice 
interactionprocesses have been extensively studied in the 
past 0000E1E1) i.e., the modulation interaction 
and the Raman interaction. The modulation interaction 
stands for the direct modulation of the inter-atomic spin- 
spin interaction due to lattice vibrations. The Raman in- 
teraction represents the second-order transition between 
the Kramers degenerate ground states of an ion through 
the intermediary of the excited states caused by the time- 
dependent variation of the crystalline field. We note 
that, for Tb 3 Ga 5 0i2, magnetic ordered states do not oc- 
cur down to very low temperature Q.2K [J^, suggest- 



ing that the inter-atomic spin-spin interaction is possibly 
very weak. We will neglect the spin-spin interaction, and 
focus on the Raman interaction. The Raman interaction 
is known to dominate the spin-lattice relaxation in many 
ionic insulators [101111. 

The electron spin-orbit coupling of a rare-earth ion is 
usually stronger than the crystalline field. By consid- 
ering the spin-orbit coupling and intra-atomic Coulomb 
interaction first, the spins and orbital angular momenta 
of the outer-shell electrons shall form a total angular mo- 
mentum J. The ground states are a J multiplet, which 
further split in the presence of the crystalline field • 
In order to develop a transparent theory, we study a rel- 
atively simple case that was often considered in litera- 
tures 0011 El El We assume that all the ground- 
state degeneracies of the ion except for the Kramers one 
are lifted by the crystalline field in such a manner that 
the energy difference E\ between the lowest excited states 
and the ground states is greater than the Debye energy 
Huj d . In this case, one can obtain for the Raman inter- 
action ^1 V = K o J2 m s m ■ Sl m , where K is a pos- 
itive coupling constant, the 1/2 isospin s m describes the 
ground-state Kramers doublet, and ft m = u m x p m is the 
center-of-mass angular momentum of the unit cell with 
u m and p m the center-of-mass displacement and momen- 
tum. The Raman interaction is taken to be isotropic 
by virtue of the cubic lattice symmetry. At low tem- 
peratures, where thermal excitation of the ion into the 
excited states is virtually impossible, the primary contri- 
bution to the magnetization comes from the ground-state 
doublet. The Kramers doublet carry opposite magnetic 
moments, which split in the presence of a magnetic field 
B and give rise to a magnetization M. It is easy to 
prove that the ensemble average of the isospin is propor- 
tional to the magnetization, namely, (s m ) = cM with c 
the proportionality coefficient . Under the mean- field 
approximation, the Raman interaction reduces to 



V = K ^2 M ■ SI, 



(2) 



where K — cKo, and KM has the units of frequency. 
While Eq. (0 is obtained in the relatively simple case, 
it represents a minimal form of possible interaction be- 
tween magnetization and phonons, and may also serve as 
a reasonable hypothesis for a phenomenological theory of 
the PHE in general paramagnetic dielectrics, similar to 
the theory of the AHE 0. 

In order to illustrate the basic mechanism of the PHE, 
we consider the vibration of the m-th unit cell in the 
phonon modes of a given wave vector q, as shown in Fig. 
1. The sound speed c T of the TA phonons is generally 
smaller than that of the LA phonons c L [l4|, meaning 
that the vibration of the unit cell in the transverse di- 
rections is lower in energy than that in the longitudinal 
direction. This yields a correlation between the angular 
momentum fl m and phonon momentum Tiq that parallel 




FIG. 1: Lowest-energy configuration between the phonon 
wave vector q and the center-of-mass angular momentum fl m 
of a unit cell. Here, R m represents the equilibrium posi- 
tion of the unit cell, which is taken to be at the origin. Top 
right panel shows the dispersion relations of the LA and TA 
phonons in the Debye model. Inset at the bottom right is a 
hypothetical level graph of a paramagnetic ion with |±) as 
the ground-state doublet. 



and antiparallel alignments between J7 m and Tiq are en- 
ergetically favorable. We may regard J7 m as an internal 
"spin" degree of freedom of the phonons. Then the corre- 
lation between i~2 m and Tiq plays the role of a "spin-orbit" 
coupling, effectively similar to the Luttinger spin-orbit 
coupling 0. Besides, the Raman interaction Eq. is 
a counterpart to the exchange coupling between carrier 
spins and magnetization in the AHE systems 1]. With 
these two essential ingredients, we can expect the PHE to 
occur in the present system, as a bosonic analogue to the 
AHE 0. In Fig. 1, the lowest-energy configuration be- 
tween fl m and Tiq is shown. It has been considered that, 
when M ^ and the Raman interaction arises as a weak 
perturbation, the degenerate TA modes split into two 
circularly polarized modes with nonvanishing i~2 m paral- 
lel and antiparallel to Tiq, respectively, and the one with 
Sl z m < is lower in energy. The splitting occurs in this 
way because |f2 m | is maximized in the circularly polar- 
ized modes, which is favorable for minimizing the Raman 
interaction. 

We start the rigorous calculation by rewriting the Ra- 
man interaction Eq. in the second-quantization rep- 
resentation 

* , V ^QO" 

C|,(7, (T * ^ 

(3) 

where A qaa i = —ihKM. - (e qa x e qa i ) with e qa the phonon 
polarization vector, and the phase convention e*_ qa = e qa 
is adopted. The thermal Hall conductivity K xy can be 
calculated by using the Kubo formula [l6| 



1 J y 



h/k B T r oo 



d\ dt(j*(-i\)jy(t)) , (4) 



where is the a;— component of the energy flux oper- 
ator J E of the phonons, and J E (i) = e iHt / h J B e- im / A . 
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The complete ex pres sion for J E for free phonons was de- 
rived by Hardy |17| . which consists of some harmonic 
terms and cubic terms. In the same essence as the har- 
monic approximation to the real lattice Hamiltonian in 
the fundamental phonon theory 5], it is sufficient to re- 
tain the harmonic terms at low temperatures. We extend 
the derivation 

El 

to include the contribution from the 

Raman interaction 
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(R m -R„)$ Q/3 (R m -R.n)< 



(•5) 



where and v" n with a = x, y and z are the a- 
components of the center-of-mass displacement u m and 
velocity v m of the m-th unit cell, respectively, and 
<j> a/3 (R m — R„) are the stiffness matrix elements of the 
lattice with R m the equilibrium position of the unit cell. 
The center-of-mass velocity is v m = p m /M c +K(M.xu m ) 
with M c the mass of a unit cell, where the second term 
originates from the Raman interaction. By using the ba- 
sic relation $^(q) = M c £ CT (w 2 ff e" e!£), the energy 



flux is derived to be J F 
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jl 1 ' , where 
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x (a_ 9ff + a\ a ){a q(T > + a^ q<yl ) , 



j q<j<j' = fWqaSvatV ' q U) qa + - (uj q(7 ~ UJ qa ,) 

x ^[(V g e-)e^,-e-(V ? e^)] 



(6) 



(7) 



(8) 



Here, 3^ comes from the Raman interaction. 

The energy scale hKM of the Raman interaction is 
generally much smaller than the Debye energy Tioj-d as 
well as the average energy difference between the LA and 
TA phonon branches, so that a perturbation treatment 
of Eq. 10 suffices. Since the two TA branches are de- 
generate, according to the degenerate perturbation the- 
ory, we need to find suitable phonon polarization vec- 
tors, for which the direct coupling between the two TA 
branches A 9 i2 = A* 21 vanishes. Such polarization vec- 
tors are obtained as e g o = «q, e g i = (9 q + i(p q )/^/2 and 

= (&q — i<£ q )/V2, where q, 8 q and (p q are the unit 
vectors associated with wavevector q = (q, 9 q , (p q ) in the 
spherical polar coordinate system. We notice that e q i 
(e q 2) for the TA phonons is a superposition of two orthog- 
onal linear polarization vectors 6 q and (p q with a fixed 
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FIG. 2: Thermal Hall conductivity function of 

KM/uJo for two different values of 8 = c L /c T . Inset: k xv 
as a function of 8 for KM = 0.001lj d . Here, T = 5.45K, 
hu n = 0.05eV and u = (12A) 3 . 



phase shift 90° (—90°) between them. This indicates that 
the TA phonons indeed split into two circularly polarized 
branches upon the perturbation of the spin-lattice inter- 
action. The standard nondegenerate perturbation theory 
is then applied to calculate the many-particle eigenstates 
of the system to the linear order in the Raman interac- 
tion. By substitution of the eigenstates into Eq. 10} , we 
consequently obtain for the thermal Hall conductivity 



-fk B KM 
2^ 2 c, 



h 



6d/T 



-dx 



1 



(9) 



where 7 = (5 - S)(l + <5) 4 /[4<5 2 (9 + 188 3 )^ 3 } with 8 = 
c L /c T , c s is the average sound speed defined by 3/c 3 = 
(1/c 3 + 2/4) and 6 D = huj /k B = (67r 2 /i/ ) 1 / 3 ^c s /fc B is 
the Debye temperature with i/q the volume of a unit cell. 
The next-order correction to n xy will be of the order of 
(KM/uj d ) 3 , which can be expected to be extremely small. 

In Fig. 2, the calculated n xy at T = 5.45K is plotted as 
a function of KM/uj d for two different values of 8. Here, 
we set Tiuj u = 0.05eV and v = (12A) 3 , which correspond 
to D ~ 580K and c s ~ 850m/s. These are either known 
values for Tb 3 Ga 5 0i2 or typical values for rare-earth gar- 
nets with similar structures ^3 ■ ^ is clear that the ther- 
mal Hall conductivity n xy is linear in the magnetization 
M. When the applied magnetic field B is relatively weak, 
M varies linearly with B, and so does K xy . We can ex- 
pect that, while K xy remains linear as a function of M, 
it may possibly become nonlinear as a function of B at 
very strong magnetic field. It is worthwhile to test this 
prediction in experiment. In the inset of Fig. 2, n xy as a 
function of 8 is plotted. K xy decreases with increasing 8. 
From Eq. ©, we see that n xy will change sign at 8 = 5. 
We notice that Eqs. © and Q make comparable oppo- 
site contributions to n xy . Their competition accounts for 
the nontrivial dependence of K xy on 8. In most materi- 
als, the typical values of 8 are around 2 where K xy 
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FIG. 3: K xy as a function of normalized temperature for 
KM — O.OOIljd and 8 = 2 (solid line). Other parameters 
are chosen to be the same as in Fig. 2. The dashed line is a 
linear fit to n xy at low temperatures. 



is always positive. 

In Fig. 3, k xy for a fixed value of KM is shown as 
a function of normalized temperature T/0 D . At low 
temperatures T <C D , we can obtain from Eq. @ 
n xy ~ (1.6'i'jk B KM/2Tr 2 c s )(k B T/h), where the numeri- 
cal identity J2 n >i n ~ 2 — 1-64 nas been used. From Fig. 
3, we see that the above linear dependence of n xy on tem- 
perature is well valid for T < O.20 D . From experimental 
point of view, it is relatively convenient to consider fixed 
magnetic field instead of fixed magnetization. For para- 
magnetic materials and relatively weak magnetic field, we 
may use the Curie's law M cx B/T. As a result, at fixed 
magnetic field, n xy approaches a constant at low temper- 
atures. At relatively high temperatures T > O.20 D , it 
decreases with increasing temperature. 

We can show that the present theory reproduces the 
essential characteristics of the experimental data In 
the experiment Q, a constant longitudinal heat current 
is driven through the sample, and the temperature dif- 
ference AT in the transverse direction across the sam- 
ple is measured. Firstly, the present theory predicts 
Kxy > for typical values of S, meaning that, if the 
magnetic field is along the z-axis and a heat current 
is driven along the £-axis, the Hall heat current will 
flow along the y-axis, which is consistent with the ex- 
perimental relation among the sign of AT and the di- 
rections of the magnetic field and the driving heat cur- 
rent Q. Secondly, AT is observed to be linear in the 
magnetic field B |4J. By analogy with the charge Hall 
effect, AT is proportional to the thermal Hall resistivity 
r xy = k xv /(k 2 xx + n 2 xy ) ~ K xy /i4 x , where the longitudi- 
nal thermal conductivity k xx depends weakly on B |4J . It 
follows that the experimental thermal Hall conductivity 
K xy oc B. In the present theory, K xy oc B at relatively 
weak magnetic field, which agrees with the experimen- 
tal observation. Thirdly, if KM is taken to be about 
l(T 4 w D , we see from Fig. 2 that n xy is lCT 5 -l(r 4 WK/m, 
which will be comparable to the experimental value 
4.5 x lCT 5 WK/m at B = IT and T = 5.45K. Here, 



the experimental value of K xy is deduced from the re- 
ported thermal conductivity lO^WK/m and 
Hall angle k xv /k xx = 1 x 10~ 4 rad T _1 0. To estimate 
the Raman coupling strength, we can use the relation 
KM — i£o|(s m )|. For the experimental material, while 
the average of the isospin is affected by many factors, 
such as the actual crystalline field and multiple para- 
magnetic ions in a unit cell, we may conservatively as- 
sume that the effective |(s m )| is of the order of O.lh at 
B = IT and T = 5.45K, where the Zeeman energy is 
already comparable to k B T . Thus Kq is estimated about 
Kq ~ 10~ 3 LL} o /h or h 2 Ko ~ 1cm -1 , which can be found 
to be within the possible range of the Raman coupling 
strength in ionic insulators US- 

In summary, we have developed a simple theory based 
upon the Raman spin-lattice interaction for the PHE, 
which has been shown in good agreement with the experi- 
mental data. Further investigation by including disorder 
effect is highly desirable. Since the PHE is similar to 
the AHE, we expect that the thermal Hall conductivity 
might be insensitive to weak disorder scattering 0. 
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